N91-20510 

MODELING OF CONNECTIONS BETWEEN SUBSTRUCTURES 

Thomas G. Butler 
BUTLER ANALYSES 

The focus of this paper is on joints that are only par- 
tially connected such as slip joints in bridges and in ship 
superstructures or sliding of a grooved structure onto the rails 
of a mating structure as shown in the sketch. 



In substructure analysis it is desireable to organize 
each substructure so as to be self contained for purposes of 
validity checking. If part of the check is to embrace a connec- 
tion, then all of the elements of the interface that it sees in 
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it 3 mate should be included within its model. In the case of the 
aroove/ rail structure, shown above, it will enhance the checking 
if the rail points, to which the shoe points will connect, are 
duplicated in the substructure with the shoe. Thus a complete 
job of checking out the shoe substrucure can be done in Phase 1 
with statics and eigenvalues and not protract the checking proce- 
dure of basic substructures into Phase 2. 

To implement such a scheme, referring to the sketch, 
points R1 & R2 are included in the shoe model. The connection 
from S3 to R1 and from S4 to R2 are made in Phase 1 and now 
become available for complete chekout of the shoe substructure, 
including its mating with the rail. To make this example 
general, postulate that the planes through the four points are 
not parallel to the coordinate planes. In eriect there are 
offsets. Generally, one likes to plan to avoid having 
out-of -plane offsets, butexigencies do crop up which forces the 
analyst to face up to such realities. Often such interface 
connections involve MFCs or elastic ties. In any case a 
requirement of Substructure Analysis is that points that are to 
be connected in Phase 2 must be available in Phase 2; i.e. they 
cannot be condensed out or constrained out in Phase 1. 
Therefore, if an MPC is used, the connecting points must be the 
independent degrees of freedom in the MPC relationship. 

The needs of this joint are that there will be no rela- 
tive translation in either transverse direction and no relative 
rotation about the long axis of the rail. In terms of the 
indicated coordinate system, translations in x and z directions 
must be constrained together and rotations about y must be 
constrained together. Just a pair of connecting points will be 
used herein to carry on the discussion. A sketch will be used to 
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assist in the discussion of making' the connection by means of 
multi-point constraints. 



Include rail point 243 in the Shoe Model. When Phase 2 COMBINE 
operation is invoked, NASTRAN will recognize that rail 243 = shoe 
243. As remarked above, since point 243 is going to be commanded 
to connect in Phase 2, it must therefore be an active available 
point for joining; and must therefore be an independent point in 
an MPC relationship. Now following the needs of this joint, 
constrain point 7013(X,Z,$) to 243<X,Z,$). The constraint 
equations for translations in X and Z are: 


7013 (X) = 243 ( X ) - c x 243($) + b x 243(¥) 

7013<Z) = 243 ( Z ) + a x 243($) - b x 243(0). 

But 243 (¥) and 243(0) are rail rotations which are not sensed by 
the shoe. If 243(4,6) are included in the shoe model they would 
be independent shoe rotations which will engage in the MPC 
relationship but would have no elastic path out to other parts of 
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the shoe. Thus, if nominal mass were added to these rail points 
to keep the eigenvalue matrix from being singular, an eigenvalue 
check for rigid body modes would show the shoe model to fail. 
One might argue, why not leave the rotations in until they are 
connected during COMBINE, then they are no longer disjoint. I 
cannot afford to leave the 243(4,6) rotations in the shoe model, 
because after connecting with the rail these rail rotations must 
not be transmitted back to the shoe. Moments in the shoe/ rail 
configuration about the two transverse axes are produced only by 
couples of forces not by local rotational bending. This rules 
out the use of MPC'.s during Phase 1 in this case. There are 
other cases of connections between substructures in which MPC's 
in Phase 1 would work. The case in which there were no 
transverse offsets would work. A NASTRAN run of a simple model 
demonstrates these results in Appendix A. 

The alternative is to make a stiff elastic connection, 
but not so stiff as to cause matrix ill-conditioning. If a bar 
instead of elastic scalars is used, it will be modeled so as to 
be fully connected in all 6 degrees of freedom at the shoe end, 
but only partially connected at the rail end. At the rail end it 
must allow for sliding along the rail and not transmit rotations 
to the shoe about the rail transverse axes. This implies that 
pin flags must be used at the rail and to inhibit these freedoms. 

This stiff bar connection can be implemented the wrong 
way or the right way. One gets trapped into modeling the wrong 
way by forgetting that pm flags are applied to bar coordinates 
not to the displacement coordinates. I fell into this trap and 
will show you what happens. Then I will follow it up with the 
correct way to model it. 
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BAR CONNECTION 
WRONG WAY 



Include the rail grid points in the shoe model and apply SPC's at 
the GRID level in d.o.f.s 2,4,6. Connect the shoe point to the 
rail point with a stiff bar. Note that the connection from shoe 
GP to rail GP produces bar coordinates that are skewed with 
respect to the displacement coordinates. Thus when bar element 
coordinate 2 is pinned, a component of force still develops at 
the rail end ot the bar in the Y displacement coordinate direc- 
tion, and so the eigenvalue check for rigid body modes fails once 
again. The listing in Appendix B of a simple model, incorporat- 
ing this wrong approach, shows the constraint forces in the rigid 
body modes in freedoms Tl, R1 , & R3 to be non-negligible . Then 
the elastic mode shows large constraint forces in these freedoms. 
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BAR CONNECTION 
RIGHT WAY 



Offset 

A 


Offset the bar at the shoe end so as to terminate the bar at the 
rail end so as to be perpendicular to all displacement coor- 
dinates at the rail end. This connection passes the eigenvalue 
check for rigid body modes. Appendix C is a listing of a simple 
demonstration problem of the joint modeled the fight way. Note 
that the constraint forces in freedoms T2 , R1 , & R3 are negli 
gible in rigid body modes as well as in elastic modes. 
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CONCLUSIONS 

This paper has demonstrated that complete checkout of a 
basic substructure can be done under the special circumstances of 
a sliding connection with offsets. Stiff bar connections make 
this possible so long as the bar coordinates are aligned with the 
displacement coordinates at the sliding surface. 
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APPENDIX A 

RUN WITH MPC CONNECTION 
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FREE-BODY MODAL STUDY OF SHOE/RAIL CONNECTIONS WITH MPC 
ALIGNED WITH OFFSETS. RAIL END INDEPENDENT. 
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MODELING OF CONNECTIONS BETWEEN SUBSTRUCTURES 


APPENDIX B 

RUN WITH WRONG BAR CONNECTION 
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GRID 
POINT 
ID. TYPE 

71 G 

72 G 


POINT 

SINGULARITY 

ORDER 

1 

1 


SINGULARITY TABLE SPC 
LIST OF COORDINATE COMBINATIONS THAT WILL 
STRONGEST COMBINATION WEAKER COMBINATION 
5 4 

5 4 


0 MPC 0 
REMOVE SINGULARITY 
WEAKEST COMBINATIO 


6 ROOTS BELOW 1.973921E+01 


EIGENVALUE 


ANALYSIS 


SUMMARY 


(INVERSE POWER METHOD) 


NUMBER OF EIGENVALUES EXTRACTED 

NUMBER OF STARTING POINTS USED 

NUMBER OF STARTING POINT MOVES 

NUMBER OF TRIANGULAR DECOMPOSITIONS . . . 

TOTAL NUMBER OF VECTOR ITERATIONS .... 

REASON FOR TERMINATION 

LARGEST OFF-DIAGONAL MODAL MASS TERM . . . . 


MODE PAIR 

NUMBER OF OFF-DIAGONAL MODAL MASS 

TERMS FAILING CRITERION 

(* 1 OR MORE ROOT OUTSIDE FR. RANGE. 
SEE NASTRAN U.M. SECTION 2.3.3) 


REAL EIGENVALUES 


MODE EIGENVALUE CYCLIC 
NO . FREQUENCY 


GENERALIZED 

MASS 


GENERALIZED 

STIFFNESS 


A 


AAAAAAAAAA ****************************************** 


NASTRAN INFORMATION MESSAGE 3308, LOWEST EIGENVALUE FOUND * 

* AS INDICATED BY THE STURM'S SEQUENCE OF THE DYNAMIC MATRIX * 

a (THIS MESSAGE CAN BE SUPPRESSED BY DIAG 37) * 


1 -7 . 963921E-08 

2 2 . 852769E-08 

3 9 . 751177E-08 

4 1 . 989545E-07 

5 2 . 305249E-07 

6 3 . 321909E-07 

7 7 . 002442E+06 


4 . 491420E-05 
2 . 688150E-05 
4 . 969911E-05 
7 . 098997E-05 
7 . 641507E-05 
9. 173054E-05 
4 . 211578E+02 


2 . 016630E-01 
1 . 250841E-01 
1 . 946527E-01 
1 . 088651E-01 
2. 362906E-01 
1. 223466E-01 
1 . 515952E-01 


A A 

-1 . 606029E-08 
3 . 568360E-09 
1 . 898093E-08 
2. 165920E-08 
5 . 447087E-08 
4. 064244E-08 
1 . 061536E+06 


35 


ID OFFSET, CONNECT 

APP DISP 

SOL 3 , 0 

DIAG 8,21,22 

TIME 10 

CEND 

FREE-BODY MODAL STUDY OF SHOE/RAIL CONNECTIONS WITH BARS JAN 20,1991 PAGE 2 
WRONG HAY WITH CONNECTOR BAR SKEWED TO RAIL. 

CASE CONTROL DECK ECHO 

TITLE = FREE-BODY MODAL STUDY OF SHOE/RAIL CONNECTIONS WITH BARS 

SUBTITLE = WRONG WAY WITH CONNECTOR BAR SKEWED TO RAIL. WANT 3 RB MODES. 

OUTPUT 

DISP = ALL 

MPCFORCES = ALL 

ELFORCES = ALL 

SPCFORCES = ALL 

SUBCASE 1 

LABEL = BARS PINNED AT RAIL END. NO OFFSETS AT SHOE END. 

METHOD = 3 
BEGIN BULK 

SORTED BULK DATA ECHO 


— 1 

++-4- 2 + + +■ 

3 

+++4+++ 

5 

+++ 6 +++ 

7 

+++8-M-+ 9 

+++ 10 +++ 

CBAR 

1 

1 

71 

7 ° 

1.0 

1.0 

0.0 

SHOE 

CBAR 

2 

T 

4 

71 

13 

14 



+TIE UP 

+TIE UP 


246 







CBAR 

3 

2 

72 

14 

13 



+TIE DWN 

+TIE DWN 

246 







CMASS2 

121 

0.1 

13 

4 




31THETA 

CMASS2 

141 

0.1 

14 

4 




14THETA 

CMASS2 

711 

0.1 

71 

4 




71THETA 

CMASS2 

712 

0.1 

71 

5 




71 PHI 

CMASS2 

721 

0.1 

72 

4 




72THETA 

CMASS2 

722 

0.1 

72 

5 




72PHI 

EIGR 

3 

INV 

0.0 

1.0 

6 

6 

3 1.-3 

+ALLMODE 

+ALLMODEMAX 








GRID 

13 

0 

3.0 

2.0 

2.0 

0 

246 

RAIL1PT 

GRID 

14 

0 

3.0 

17.0 

2.0 

0 

246 

RAIL1PT 

GRID 

71 

0 

0.0 

0.0 

0.0 

0 


SHOE1PT 

GRID 

72 

0 

0.0 

15.0 

0.0 

0 


SHOE1PT 

MAT1 

1 

3. +7 


0.28 

2.4-4 




PARAM 

COUPMASS7 







PBAR 

1 

1 

1.0 

1.0 

1.0 

1.0 



PBAR 

ENDDATA 

2 

1 

100.0 

100.0 

100.0 

100.0 
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1 


SUBCASE 1 

REAL 


EIGENVALUE = -7 . 963921E-08 
EIGENVECTOR NO. 


PT ID. T1 T2 

13 3 . 308091E-01 0.0 

14 -6. 346744E-01 0.0 

71 3.910763E-01 -5 . 238029E-01 


T3 . R1 R2 R3 

5 615981E-01 0.0 9.189782E-02 0.0 

-5 552183E-01 0.0 9.189782E-02 0.0 

1 ! OOOOOOE+OO -1 . 011211E-01 6.503133E-02 6.436E-, 


N 0 


iIJSSSESi -i.ximis-oi 6.503133.-02 6 . 436.-2 

EIGENVALUE = 2.852769E-08 REAL EIGENVECTOR 

T3 R1 R2 M 

o 554041E-01 0.0 - 1 . 205854E-01 0.0 

- SHEli J tiii h\sz 


PT ID. T1 
1 -1 . 806720E-01 

14 -7 . 829153E-01 


T2 

0.0 

0.0 


EIGENVALUE = 9.751177E-08 REAL EIGENVECTOR NO 


PT ID. T1 

13 -7 . 289532E-01 

14 3 . 520157E-01 

71 -8.096893E-02 

72 


T2 T3 

0 0 2 . 853720E-01 

0.0 2.345846E-01 

-6. 163144E-01 -7 . 029006E-02 


R1 

0.0 

0.0 

- 3 . 385834E-03 


lioOOOOOE+OO -6. 163144E-01 - 1 . 210776E-01 - 3 . 385834E-03 


R2 

-1.817658E-01 
-1.817658E-01 
- 2 . 055022E-01 -7.2E- 
- 2 . 055022E-01 -7.2E- 


R3 

0.0 

0.0 


EIGENVALUE = 1.989545E-07 REAL EIGENVECTOR 


N 0 


PT ID. T1 

13 3 . 177751E-01 

14 -4.354112E-01 

71 1 . OOOOOOE+OO 

72 


'T'o T3 RI 

o 0 1.764639E-02 0.0 

0*0 3 . 301453E-01 0.0 

-9 . 989289E-02 -9 . 057981E-01 2.083326E-02 


r: 


R3 


■3. 130602E-01 0.0 


-3 . 180602E-01 
. 929949E-01 


0.0 

5. 02E-2 


i:SSS*?£oS *•«»-* 


EIGENVALUE = 2.305249E-07 


REAL EIGENVECTOR 


N 0 


PT ID. T1 T2 T3 

13 9.812109E-01 0.0 1. 000000E+00 

14 8 . 707730E-01 0.0 5.573768E-01 

71 7 983168E-01 4.514050E-01 8.229362E-01 


R1 R2 R3 

0 o -1 . 272332E-02 0.0 

0*0 -1 . 272332E-02 0.0 

2.950821E-02 - 2 . 407710E-02 7.363E-3 


N 0 


71 7 . 983168E-01 4.514050E-01 -^^lE-oi -2 407710E-02 7.363E-3 

72 6.878789E-01 4.514050E-01 3.803130E-01 -2.950821E 02 i . 

EIGENVALUE = 3.321909E-07 REAL EIGENVECTOR 

mn R1 R2 

PT ID. T1 T2 nnnnfiOF+OO 0 0 - 3 . 662173E-01 0.0 

13 -7.600989E-01 0.0 3*?S5??2E-01 00 -3 . 662173E-01 0.0 

I £3253 «SW :I:iIlItlI:Si 


EIGENVALUE = 7.002442E+06 

PT ID. T1 T2 

13 -6.670417E-01 0.0 

14 

71 

72 


REAL EIGENVECTOR 


N 0 


5.339078E-01 0.0 
6 . 416196E-0 1 1 . 446323E-03 
-5 . 595077E-01 7.309780E-03 


T3 R1 R2 R3 

1 . OOOOOOE+OO 0.0 ■ 6,548 °?!r’!l n'n 

-R 777O38E-01 0.0 5.701413E010.) 

-9 ! 644380E-0 1 2.011336E-01 '5 • 263 158E-01 J-^E-l 
8.305104E-01 5.448166n-0^ o.091^.5E 01 4.b^ 
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CONSTRAINT 


EIGENVALUE = -7 . 963921E-08 FORCES OF SINGLE-POINT 

PT ID. T1 T2 T3 R1 R2 

13 0.0 2 . 806213E+01 0.0 -2 . 590351E+01 0.0 

14 0.0 -1 . 530662E+01 0.0 -4 . 709728E+00 0.0 


EIGENVALUE 

= 2 . 352769E-08 

PT ID. T1 

13 0.0 

14 0.0 

T2 

2 . 551 103E+00 
-2. 551103E+00 

EIGENVALUE 

= 9 . 751177E-08 

PT ID. T1 

13 0.0 

14 0.0 

T2 

-1 . 275551E+00 
5. 102206E+00 

EIGENVALUE 

= 1 . 989545E-07 

PT ID. T1 

13 0.0 

14 0.0 

T2 

9. 566635E-01 
-5. 102206E+00 

EIGENVALUE 

= 2 . 305249E-07 

PT ID. T1 

13 0.0 

14 0.0 

T2 

1 . 020441E+01 
2. 793963E-09 

EIGENVALUE 

= 3 . 321909E-07 

PT ID. T1 

13 0.0 

14 0.0 

T2 

-2 . 551103E+00 
-2. 551103E+00 

EIGENVALUE 

= 7 . 002442E+06 

PT ID. T1 

13 0.0 

14 0.0 

T2 

-3 . 846585E+01 
-1 . 149591E+02 


FORCES OF SINGLE-POINT 

T3 R1 R2 

0.0 1 . 177432E+00 0.0 

0.0 - 1 . 118560E+01 0.0 

FORCES OF SINGLE-POINT 

T3 R1 R2 

0.0 1 . 876532E+00 0.0 

0.0 1 . 103843E-01 0.0 

FORCES OF SINGLE-POINT 

T3 R1 R2 

0.0 3 . 017170E+00 0.0 2 

0.0 4 . 121012E+00 0.0 2 

FORCES OF SINGLE-POINT 

T3 R1 R2 

0.0 1.964840E+01 0.0 

0.0 -5 . 077676E+00 0.0 

FORCES OF SINGLE-POINT 

T3 R1 R2 

0.0 -3 . 532296E+01 0.0 

0.0 -1 . 236304E+01 0.0 

FORCES OF SINGLE-POINT 

T3 R1 R2 

0.0 -2 . 479790E+04 0.0 

0.0 2 . 221343E+04 0.0 


R3 

-1. 726900E+01 
-3 . 139819E+00 

CONSTRAINT 

R3 

7 . 849547E-01 
-7 . 457070E+00 

CONSTRAINT 

R3 

1. 251022E+00 
7. 358950E-02 

CONSTRAINT 

R3 

;. 011446E+00 
!. 747341E+00 

CONSTRAINT 

R3 

1 . 309893E+01 
-3 . 385117E+00 

CONSTRAINT 

R3 

-2. 354864E+01 
-3 . 242024E+00 

CONSTRAINT 

R3 

-1 . 653193E+04 
1 . 480896E+04 
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MODELING OF CONNECTIONS BETWEEN SUBSTRUCTURES 


APPENDIX C 

RUN WITH RIGHT BAR CONNECTION 
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